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Abstract 

<N ■ 

O ■ We present the pion-pion elastic scattering amplitude at finite temperature to 

one-loop in Chiral Perturbation Theory. The thermal scattering amplitude properly 
defined allows to generalize the perturbative unitarity relation to the T / case. 
Our result provides a model independent prediction of an enhanced pion-pion low- 
energy phase shift with the temperature and it has physical applications within the 

t^J- | context of Relativistic Heavy Ion Collisions. 

^ ! PACS: ll.10.Wx, 12.39.Fe, 11.30.Rd, 25.75.-q. 

> : 

2 : 1 Finite Temperature Pion Scattering 

m : 

There is a growing interest in the Quark-Gluon-Plasma produced in Relativistic Heavy 
Ion Collisions (RHIC) and its evolution into hadronic matter. After the chiral phase tran- 
sition, this plasma cools down rapidly, mainly into a pion gas, which requires a description 
Ph| of QCD at temperatures below the chiral phase transition. With this aim we turn to Chi- 

ral Perturbation Theory (ChPT) [1, 2] which is the low-energy effective theory of QCD. It 
follows from the identification of pions as the Nambu-Goldstone bosons associated to the 
QCD SU (2) l x SU {2)r — > SU (2) l+r spontaneous chiral symmetry breaking and has been 
already applied to describe static properties of the low-T pion gas [3]. Within ChPT, we 
will analyze the thermal tcti scattering, which is relevant for the gas description in terms 
of transport equations, for the evolution of the expanding fireball or for estimating the 
freeze out time (when the scattering mean time is larger than the collective expansion 
time scale). Pion scattering in the a channel has also been related to chiral symmetry 
restoration [4, 5]. The thermal amplitude in the vector channel has been studied within 
a modified Nambu-Jona Lasinio (NJL) model [6] explicitly including the prm vertex and 
analyzing the complications due to a non-physical constituent quark-antiquark threshold 
near the p mass. Within ChPT only the scattering lengths are available at finite T [7] (for 
a NJL calculation see [8]). In this work we calculate the whole nir elastic scattering ampli- 
tude at finite temperature to one loop in ChPT, thus valid at energies and temperatures 
well below A x ~ 1.2 GeV, both considered as 0(p) in the chiral expansion. 

A further motivation for our analysis is one of the most striking observations of RHIC: 
the modification of the dilepton spectrum with respect to the vacuum hadronic emission 
[9]. This is usually interpreted as an in-medium modification of the p mass and width 
[10, 11]. Here we will obtain a completely model independent and systematic description 
of the low energy tail of this resonance, in terms of tttt scattering. Let us remark that 
we are taking into account the details of the pion interactions, not only their distribution 
function [11] whose effect comes from the imaginary part of the amplitude, as we will see 
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below. The description of the full resonance is beyond ChPT. Nevertheless, we show here 
that our amplitude satisfies a thermal generalization of perturbative unitarity, which is 
the starting point [12] of the successful nonperturbative unitarization methods used to 
study resonances like the p and a at T = [13]. 

2 The T > Pion Scattering Amplitude 

2.1 Definition and interpretation of the thermal amplitude 

Let us start by discussing briefly what we mean by thermal scattering amplitude in what 
follows and why it contains interesting physical information. For our purposes, it will be 
enough to assume that the system is in thermal equilibrium and to neglect finite baryon 
density effects (as in an ideal RHIC in the central rapidity region). 

In Thermal Field Theory, the word "amplitude" is often used to denote a thermal 
Green function rather than an S'-matrix element, whose definition in a thermal bath 
is rather subtle. At T = the S-matrix element is related to a time-ordered product 
through the LSZ reduction formula. Such a time-ordered product can be calculated at 
T / 0, where the time components are ordered along a complex contour C [14]. In 
the imaginary-time formalism (ITF) C = [0, —i/3] so that the time arguments become 
Euclidean r = it e C, with (3 = T _1 and t being the Minkowski time. The Feynman 
rules in the ITF are simply those at T = but replacing all the zeroth components of 
momenta by a discrete set of frequencies q° — * iuj n = 2mnT and the loop integrals by 
Matsubara sums, i.e., {2n)~ l f dq° — > %T X^L-oo- Since the vertices remain unchanged 
the set of diagrams is exactly the same as for T = 0. The ITF is appropriate to calculate 
thermodynamical quantities, but it can also be used to analyze dynamical processes with 
real energies. In fact, carrying out the Matsubara sums first, the ITF Green functions 
yield analytic functions of iuo n everywhere off the real axis [14]. Therefore, they can 
be continued analytically from discrete external complex to continuous real energies E, 
usually with the prescription iou n — > E + ie. It is worth mentioning that the thermal Green 
functions thus obtained correspond to the retarded functions in the real-time formalism 
(RTF) [15, 16] for which C includes the real axis. One could also calculate the time- 
ordered product directly in the RTF, where the external energies remain real, but the 
number of fields is doubled [14]. In this sense, let us remark that in the RTF it has 
been shown [16] that the T ^ retarded functions have a causal and analytic structure 
which allows for a description in terms of dispersion relations, unlike the T ^ RTF 
time-ordered product. Hence, the retarded functions are suitable candidates to generalize 
important properties of the T = amplitude such as unitarity, as we will see below. 

Therefore, back to our present case, consider the T = S'-matrix element correspond- 
ing to 7T7T scattering, written through the LSZ formula in terms of the Green function 
of the time-ordered product of four pion fields. The thermal scattering amplitude we 
will consider is defined with respect to such T = amplitude by keeping the pion initial 
and final asymptotic states as those for T = but i) replacing the Green function by a 
thermal expectation value in the ITF and ii) performing the usual analytical continua- 
tion into continuous real energies. As discussed above, this corresponds to the retarded 
four-point function in the RTF and, as we will see, the thermal amplitude thus defined 
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has the correct analytic structure, allowing for a suitable extension of the perturbative 
unitarity at T 7^ 0. Let us remark that the same definition for the pion scattering thermal 
amplitude has been used in [6, 7, 8]. 

An alternative point of view would have been to consider the RTF thermal reaction 
rate 1Z for the inclusive reaction 7T7r+Heat Bath — > nn + anything [17] defined as the 
thermal average of the transition probability S*S, with the T = ^-matrix element 
between an initial and final asymptotic state including the heat bath collective states. We 
remark that this is a different object than the thermal amplitude we consider here, thus 
providing information on different physical observables. Thus, while TZ can be interpreted 
as an "average of probability" when the free heat bath is taken also as an asymptotic 
state, the thermal amplitude we consider here can be viewed rather as (the retarded part 
of) an "average of amplitude" of the probability that two initial free pions scatter in the 
interacting thermal bath and two final free pions leave it. In principle, there is no obvious 
relation between the two, except for the case of two-point functions, where the imaginary 
part of the retarded RTF Green function is proportional to the difference between the 
direct decay rate, given by 1Z and the inverse one, related to TZ by detailed balance [17, 16]. 
A detailed calculation of the T ^ reaction rate for the pion scattering using the rules 
given in [17] would be very interesting but it lies beyond the scope of this work. 

2.2 The one loop ChPT calculation 

The Nf = 2 chiral lagrangian at 0(p 2 ) is that of the 0(4) nonlinear sigma model: 



where / is the T = pion decay constant in the chiral limit, U(x) = (Jl — it 2 / f 2 , n/f) is 



an 0(4) vector satisfying U T U = 1, with ir a (a = 1, 2, 3) the pion fields and \ — ( m2 > 



Here, m is the pion mass to lowest order (in the isospin limit). For simplicity, we do 
not show external sources, which are nevertheless needed to obtain the 0(p A ) relations 
between the physical m n , f n and m, f [2]. The 0{p A ) lagrangian £ 4 includes seven 
independent terms multiplied by low-energy constants, called li after renormalization. 

There is another subtlety to be borne in mind. The loss of Lorentz covariance imposed 
by the choice of the thermal bath rest frame implies that we have to separate explicitly 
the time and space-like components of any four-vector. In particular, for an elastic two- 
body collision with four-momenta k\k 2 — > k 3 k 4 , if we call S = ki + k 2 , T = k\ — k 3 
(not to be confused with the temperature) and U = ki — k±, the amplitude will depend 

— » — » — * 

differently on S , |S|, T , |T|, U and |U|. In particular, crossing symmetry on the 
Mandelstam variables s = S 2 , t = T 2 , u = U 2 no longer holds. Let us recall that at T = 
any elastic irn amplitude is related to that of 7r + 7r~ — > 7r°7r°, called A(s,t,u), by isospin 
and crossing transformations. Nevertheless, since the temperature does not modify the 
vertices, crossing symmetry still holds in terms of the S, T and U four-momenta, and we 
can still write any ixix thermal amplitude in terms of the thermal A(S, T, U; j3): 



A(S, T, U; 0) = A 2 {S, T, U) + Af (S, T, U) + Af d (S, T, U; (3) + ^(S, T, U; /3), (2) 




(1) 
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where we have followed the usual T = notation: A 2 is the £ 2 tree level contribution and 
A^° l contains the £4 tree level, depending on the fj, plus those polynomials coming from the 
renormalization of the loop integrals. Both A 2 and A v £ l are temperature independent. The 
j^uni t erm represents those contributions from loops with two propagators proportional 
to the three independent integrals Jo, J\ and J 2 defined in eq.(14) which yield the correct 
analytic structure at T = and will ensure perturbative unitarity in all channels. The 
terms collected in A 1 ^ are proportional to the tadpole integral F@ in eq.(13) and include 
all tadpole diagrams as well as terms coming from loop integrals (see Appendix A). 

Taking these remarks into account, the thermal amplitude is obtained by adding to 
the temperature independent part in [2] the following temperature dependent corrections: 



AA* a<i (S, T, U; /3) = 



f 4 { 
2(£. *)(£•*) 



-s + t + u + 



5m 2 



^-(k% + k%) f + 2 (fci - fc 2 ) • T 
I T 1 2 



> ( - + To' 



+ ^lM + (T^U) 



(3) 



AA% ni (S, T, U; (3) = 



Ij ( (. - m >) [-2^,(8) + ^ (m 2 + 3, - 2* 3 .S) - ^(sA M S)-^ A M S)) 



+ 



m 



(jfci - k 2 ) -T [tAJo(T) - 2T AJi(T)]] + m 4 AJ (T) - 2(A;° - k° 2 )m 2 AJi(T) 



-4A;°A; 2 AJ 2 (T) - (fc°fc 2 + fcife) f [tAJi(T) - 2T AJ 2 (T)] 



2(A; 1 -T)(A; 2 -T) 



(3T 2 - |f| 2 ) AJ 2 (T) - 3T tAJ!(T) + ( ^ + m 2 \f | 2 ) AJ (T) 



2fci • k 2 



[tAJ 2 (T) - *T AJi(T) + - + m 2 |T| 2 AJ (T) 



^ 4 
+ (T 



U) 



(4) 



Note that the above equations in Minkowski space-time are obtained from the Eu- 
clidean expressions by replacing k° — > — iE for all Euclidean momenta k , with E real, 
and analytically continuing the J k as showed in eq.(23). Note also that, for clarity, we 
have not simplified the ki four-momenta in terms of S, T and U. In addition, the T- 
corrections to m n and f n [3] are included in AA^ d , consistently with our definition of the 
thermal amplitude. Recall that the T = physical f n and m % are [2]: 



m 



/ 3 m 2 
32tt 2 / 5 



U = f i 



m 2 l 4 



(5) 



We have performed three different consistency checks on the thermal amplitude: First, 
we have recovered the Lorentz invariant T — > + limit of [2]. Since all the temperature 
dependence is on the tadpoles and Jk functions, we are thus checking their coefficients, 
since our calculation was carried out without Lorentz covariance. The other two checks 
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concern partial waves and the c.o.m. frame, i.e, S = or pions at rest with the thermal 
bath. Customarily, the tttt scattering amplitude is projected into partial waves ajj of 
definite isospin / and angular momentum J, and the standard definitions [2, 13] are still 
applicable at T ^ 0, but with crossing symmetry in terms of the S, T and U four- vectors. 
This procedure is equivalent to a direct calculation of the isospin amplitudes without 
using the A function, as done in [7] , whose analytic expressions for the thermal scattering 
lengths we have reobtained, thus checking the real part of our partial waves with an 
independent calculation. The third check concerns the imaginary parts of the thermal 
partial waves in the different channels. This requires a generalization of the perturbative 
unitarity relation that we present next, which, in turn, will provide a natural physical 
interpretation of the imaginary part of our thermal amplitude in terms of direct (emission) 
and inverse (absorption) processes in the thermal bath. 



2.3 Thermal unitarity 

For simplicity, we drop the 1.1 indices in the following. At T = unitarity constraints 
the partial waves, for s > Ami an d below other inelastic thresholds, to satisfy 



Imo(s) = a{s)\a{s)\ 2 , a(s) = yjl - Am 2 / s , (6) 

whereas the ChPT series only satisfies unitarity perturbatively, i.e., 

Ima 2 (s) = 0, Im a 4 (s) = a(s) \a 2 (s)\ 2 , .... (7) 

We now generalize the above perturbative relation for any one-loop calculation, includ- 
ing those within ChPT. Let us recall that for one-loop elastic amplitudes, any imaginary 
part should come from loop diagrams with two propagators and in particular from the 
one-loop integrals Jk (k = 0, 1,2) given in eq.(14). As an illustration we will analyze Jo, 
since the other cases are completely analogous. We start then by performing the frequency 
sums in the standard fashion [14]: 

T (,\n\T\- f dD ~ lq X- sis 2 l+n(s l E q )+n(s 2 E q _ Q ) 
J ° ( "' lQl T) ~ ~ J WY^ u- Sl E q -s 2 E q _ Q (8) 

where E\ = \k\ 2 +m 2 , n(x) = [eP x — l) is the Bose-Einstein distribution function and 
io = iQ = lnkTi with k integer, is the external frequency. Note that eq.(8) provides 
the analytic continuation of Jq(ou) for u> off the real axis. Taking now Im Jq(E + ie) with 
E real (see our previous discussion on the thermal amplitude) the two s\ = —s 2 terms 
cancel each other. The two cases s 1 = s 2 = ±1 require E ± E q ± E q _Q = which imply 



±E — y\Q\ 2 + 4m 2 > 0, respectively. Summarizing: 

lmJ (E + ie,\Q\;T) = tt sgn(E)6(s - 4m 2 ) (9) 

x / Wf^kr Q 11 + " w + " (Vo)l 4 m " E ' " E '- Q > 

where s = E 2 — \Q\ 2 and E real. Hence Jo(s) has a cut on the real axis starting at s = 4m 2 , 
as in the T = case. Our analysis is similar to that of a decaying particle in medium [18], 
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where the integral in eq.(8) is the thermal self-energy. In fact, the interpretation of eq.(9) 
follows from the observation that l+n(E 1 )+n(E 2 ) = [l+n(E 1 )][l+n(E 2 )\ -n(E 1 )n{E 2 ). 
Therefore, in our thermal amplitude and for positive energy, the in-medium pions increase 
the available phase space by (1+ni) (1+722) (corresponding to the thermal enhancement of 
the two pion production) and decrease it by n\n 2 (when the two initial pions collide with 
thermal pions and are absorbed) analogously to the stimulated emission and absorption 
processes discussed in [18]. Alternatively, we could have calculated in the RTF, arriving to 
the same result for the retarded self-energy, unlike the RTF time-ordered product, which 
yields the sum of the emission and absorption terms rather than their difference [16]. 

Partial waves are defined in the c.o.m. frame, where we start studying the S'-channel 
diagrams, setting Q = in eq.(9). The generalization to one-loop ChPT is straightfor- 
ward. Since C 2 contains only two derivatives, one needs loop integrals with up to four 
momenta in the numerator, whose imaginary part can be related to Jo, J\ and J 2 in 
eq.(14) (see Appendix A). The J ii2 analysis is similar to that for J , and we find 

lmJ k (E + ie,6;T) = ^p-0(s - 4m> T (£) (|)\ k = 0,1, 2 (10) 



where a T {E) = a{E 2 ) 



1 + 



exp(/3|£|/2) - 1 



(11) 



is the thermal two-particle phase space in the c.o.m. frame [18]. Note that in (10) we 
have continued analytically the J k functions to real energies as indicated in eq.(23) . 

Concerning the T and U channels, in the c.o.m. frame and for s > 4m 2 we have to 
evaluate Jfe(0, |T|) (see Appendix A) which is real, according to eq.(9). Therefore, the 
same T = derivation of eq.(7) follows for T > by replacing a — > a T) i.e.: 



Ima 2 (s)=0, Im a 4 (s; /9) = <tt(S ) |a 2 (s)| 2 , S > 2m„ 



(12) 



This generalizes the perturbative unitarity relation to T 7^ 0, providing a neat physical 
interpretation of the imaginary part of the thermal amplitude in terms of emission and 
absorption of pion pairs in the thermal bath. Note that the one-loop thermal amplitude 
considered here, in terms of retarded propagators, has the same unitarity cut as the T = 
one. In fact, it can be shown that they have the same analytic structure in the complex 
plane [12], which gives support to the idea that the thermal retarded functions provide 
the natural extension of amplitudes as far as their causal and analytic properties are 
concerned. We remark that our perturbative result (12) is consistent within ChPT, where 
the external energies and T remain small and so do the Bose-Einstein factors. For higher 
orders one could have intermediate states in the thermal bath with more than two pions, 
weighted by higher powers of n. We have explicitly checked that our final result for the 
amplitude satisfies eq.(12) for the JJ=00,11,20 channels. 

Finally, let us note that regarding the 7T7t scattering in the 11 channel as the exchange 
of a p resonance [19], the pion scattering amplitude would be proportional to the p thermal 
retarded self-energy, whose imaginary part to one loop is proportional to a T /cr and gives 
the net decay rate of the p in the thermal bath. Note that, from this point of view, the 
stimulated emission contribution, proportional to [1 + n(\E\/2)} 2 , would give the reaction 
rate discussed in [17] for the decay of a p at rest into two pions. Hence, we expect to be 
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able to extract physical information about the thermal properties of physical resonances 
arising in our nn scattering amplitude, which we have calculated without introducing 
those resonances as explicit degrees of freedom. This is indeed the shown recently 

in [12], where the thermal mass and width of the p and a are read off from the complex 
poles (in the second Riemann sheet) of the thermal nonperturbative amplitude, obtained 
from the ChPT one by unitarization. At low T, the p width grows with the a T /a factor, 
as our previous discussion suggested. This provides a relevant physical application of our 
present analysis, which is totally model independent, in contrast with the very recent 
results in [12], which maybe more direct to extract information on resonances but relies 
on the choice of unitarization method. 

3 Discussion and conclusions 

From our previous results, we have also calculated the scattering phase shifts, 5jj in 
different channels. This is a simple and easy way of representing the complex amplitude 
in terms of just one parameter, since elastic unitarity implies a = exp(£) sin(<5)/o"T- Since, 
for a given energy, the thermal phase space is fixed, the "strength of the interaction", 
i.e. the modulus of the amplitude, is given only in terms of sin S This is similar to the 
T = case, and perturbatively corresponds to St — o T (t/s) [a 2 (s) + Rea^s, (3)} . The 
results are plotted in Fig.l up to 500 MeV, roughly where ChPT is reliable. Note that the 
absolute value of all the phase shifts is increased, but each channel keeps its attractive 
or repulsive nature (the sign of the phase shifts). The phase shift increase we find is 
basically due to the Bose enhancement factor in a^, the real part and the modulus of 
the perturbative amplitude changing very slightly with T at low energies. In fact, in the 
IJ = 20 channel, |Re decreases with T for low energies (see [7]). A strong enhancement 
of | aoo 1 2 near threshold would be a signal of chiral symmetry restoration [4, 5] but we do 
not observe such effect here. We also remark that not all the thermal effects in our one- 
loop calculation can be accounted for with a thermal redefinition of /„- [5] since in that 
case, with the decreasing f n (T) given in [3], we would see a larger increase of \a\ 2 in 
all channels. Note also that our result can be alternatively interpreted as an enhanced 
effective pion interaction if we were describing them with the T = expressions, although 
bearing in mind that the main contribution comes from the phase space factor. 

The I J = 11 channel is particularly interesting, since it is dominated by the p res- 
onance, whose thermal deformation has been related to the anomalous dilepton spectra 
in RHIC [10, 11]. By analyzing this channel at low energies, we are studying the tail of 
the resonance. The low-energy temperature enhancement of the phase shift can then be 
interpreted as a stronger thermal effect of the rho, although with pure ChPT calculations 
we cannot disentangle whether it is due to a larger width or to a smaller mass. Qualita- 
tively, the chiral combination of parameters that dominate the phase shift in this channel 
is lx — l 2 ~ Tpf^/Mp, following the resonance saturation hypothesis [2, 19]. In order to 
obtain larger phase shifts at T = 0, we need somewhat larger l x — l 2 values. Thus, the 
low-energy thermal phase shifts would correspond at T = to a larger T p f 4 /Mp ratio, 
although with the T = part only we cannot mimic accurately the low- energy thermal 
behavior. Therefore, our results are consistent with most models that expect a low-T 
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Figure 1: Temperature evolution of the phase shifts 5u for I J = 11, 00, 20. 



sizable thermal increase of the rho width due to the thermal increase of phase space, with 
an almost constant thermal M p [20, 11]. Our model independent perturbative result is 
confirmed by our discussion at the end of section 2.3 and by the unitarized results in [12]. 

In conclusion, we have calculated the rnv scattering to one loop in Chiral Perturbation 
Theory at finite temperature T. Our thermal amplitude has been calculated with T = 
asymptotic states and the thermal ITF Green function, whose analytic continuation 
corresponds to the retarded RTF Green function, which has a suitable causal and analytic 
structure. Despite the loss of Lorentz covariance of the thermal formalism, any irir elastic 
amplitude can still be related to a single one by means of a generalized crossing symmetry 
in terms of the external four-momenta. We have found the thermal generalization of 
perturbative partial wave unitarity to one loop, which amounts to replacing the T = two- 
particle phase space with its thermal version. This accounts for emission and absorption 
of pion pairs inside the thermal bath, providing a physical interpretation of the imaginary 
part of the thermal amplitude and giving the expected analytic structure above threshold, 
consistently with causality and thermal unitarity. We have checked that our amplitudes 
satisfy this constraint. With these amplitudes, we have studied the low energy behavior 
of thermal 7T7T elastic scattering in a model independent formalism, finding an enhanced 
phase shift in all channels, as dictated by the thermal phase space, but with little change 
in the modulus of the partial waves. Physically, this means that one cannot mimic the 
low-T behaviour of dynamical quantities just by scaling /„. — > f n (T). 

One of the main physical observables that can be related to our thermal scattering 
amplitude are the masses and widths of the p and a mesons [12]. We have given here 
a low-T qualitative description for the case of the p. Further quantitative and detailed 
information can be obtained if these results are extended to higher energies beyond pure 
ChPT. In particular, by means of unitarization methods [13] that generate dynamically 
light resonances like the p and a mesons [12]. Hence, our results are useful as a starting 
point for the description of relevant physical processes where chiral symmetry and uni- 
tarity play a crucial role, such as the properties of those resonances in the plasma formed 
after a Relativistic Heavy Ion Collision. 
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A Loop integrals 

The following are the integrals needed for our calculation in an arbitrary reference frame 
(not necessarily the c.o.m frame, where pions are at rest with the thermal bath). The 
signature of the metric in Euclidean space is (— — ) and we will use the notation 



r d D q 1 n ^°° r d D ~ l q 



(2tt)' 



where (3 — 1/T, qo — 2irnT and Qo = 2i\kT with k integer denotes an external frequency. 
Also we denote AG(T) = G(T) — G(T = 0). First, we have the tadpole integral: 



d D q 



1 



I ft (2n) D q 2 -m 2 ' 
Then, three integrals of pion loops with an external momentum flow: 

d D q 



J k (Qo,\Q\;T) 



Qo 



, fc = 0,l,2 



h {2n) D {q 2 - m 2 ){{q - Q) 2 - m 2 
All other integrals with momenta in the numerator can be reduced to those four as 



r d D q 
J j2^W 



Qi 



Q l 



ft (2tt) d (q 2 - m 2 )((q - Q) 2 - m 2 ) 

d°q qJT 

'p (2ir) D (q 2 - m 2 )((q - Q) 2 - m 2 ) 



\Q? 
Q2 J 



QqJi + -JoQ 2 



(13) 

(14) 

(15) 
(16) 



d D q 



QoQi 



L 



ft (2ir) D (q 2 - m 2 )((q - Q) 2 - m 2 ) 
d D q qtfj 



ft (2tt) d (q 2 - m 2 )((q - Q) 2 - m 2 ) 



\Q? 

QiQjha + gijhb 



Q 2 Qo 

Q0J2 + \Ji + 



(17) 
(18) 



I 

L 



d D q 



Qoq-q 



ft (2ir) D (q 2 - m 2 )((q - Q) 2 - m 2 ) 

dPq q%q-q 

'/3 (2tt) d (q 2 - m 2 )((q - Q) 2 - m 2 ) 



= QoFp + m 2 ^ 



QiFp - m 2 -Q- \QqJi + ^ J 



Q 2 
2 



(19) 
(20) 



I 



d D q 



(q-qf 



(3 (2tt) d (q 2 - m 2 )((q - Q) 2 - m 2 ) 



{Q 2 + 2m 2 )F p + m 4 J . 



(21) 
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where 

ha{Q) 



+ 



{D-2)\Q\ 4 

(<Q 2 



((D - l)Ql + |Q| 2 ) J 2 + (D - l)Q Q 2 Ji + (^j^Q 4 + ™ 2 \Q\ 



hb{Q) 



{D-2)\Q\ 2 



+ Q${D-l) + \Q\ i \Ff> 



-Q 2 J 2 + Q Q 2 Ji 



^ + m 2 \Q\ 2 ).Jo-^F, 



The above integrals have the appropriate T — > 0, m — > limits and at T = they 
can be written just in terms of the tadpole F and J [2]. Recall that all the UV divergences 
of the Fp and J integrals are contained in the T = part, since AF and AJ k always 
contain Bose-Einstein functions vanishing exponentially for large momentum. Therefore, 
using dimensional regularization, we always separate the T = part, regularized and 
renormalized as in [2] and set D = 4 in AFp and AJ k . 

The functions Jk in eq.(14) can be analytically continued after performing the Mat- 
subara sums. All the Euclidean discrete frequencies are then continued to — iE with E 
real and continuous, i.e, 



J k (Qo,\Q\;T) - (-i) k J k (-tE,\Q\;T) 



(23) 



In the c.o.m. frame the imaginary part of the J integrals, can be calculated exactly, e.g, 
eq.(10). When the c.o.m frame |S| = T = U = 0, coincides with the thermal bath rest 
frame, only Fp, Jq(S), Jq(T), J 2 (T) have to be evaluated (note that Jo(U), ^(U) can be 
computed from the corresponding T integrals) . Specifically 



MS) = ^J (s) + f , Ji(s) = | 



J (S), Ji(T)=0 



(24) 



The real part of AJ analytically continued as in eq.(23) in the S'-channel reads 



where q(E) = V E 2 — m 2 n(E) = (e /3E — l^j and PV denotes Cauchy's Principal Value. 
Note that the above expression is valid for uo real (s = oo 2 ) and the pole of the integrand 
is only present for s > Am 2 , where Im J ^ 0. Finally, in the T(U) channels, we get 



AJ fe (|T|) 



87r 2 io q |T| 6 \2g 



k = 0,2 



where E 2 = q 2 + m 2 . The above result is valid for real |T| (or s > Am 2 for -kit scattering). 



References 

[1] S. Weinberg, PhysicaA 96, 327 (1979). 



10 



[2] J. Gasser and H. Leutwyler, Annals Phys. 158, 142 (1984); Nucl. Phys. B 250, 465 (1985). 

[3] J. Gasser and H. Leutwyler, Phys. Lett. B 184, 83 (1987); P. Gcrber and H. Leutwyler, Nucl. Phys. 
B 321, 387 (1989); A. Schcnk, Phys. Rev. D 47, 5138 (1993). 

[4] S.Chiku and T.Hatsuda, Phys. Rev. D 57, R6 (1998). 

[5] D.Jido, T.Hatsuda and T.Kunihiro, Phys. Rev. D 63, 011901(R) (2000). 

[6] Y. B. He, J. Hufncr, S. P. Klcvansky and P. Rehberg, Nucl. Phys. A 630, 719 (1998). 

[7] N. Kaiser, Phys. Rev. C 59 (1999) 2945. 

[8] E. Quack, P. Zhuang, Y. Kalinovsky, S. P. Klevansky and J. Hufner, Phys. Lett. B 348, 1 (1995). 

[9] B. Lenkeit et al. [CERES-Collaboration], Nucl. Phys. A 661, 23 (1999) 

[10] G. Q. Li, C. M. Ko and G. E. Brown, Phys. Rev. Lett. 75, 4007 (1995) 

[11] V. L. Eletsky, M. Belkacem, P. J. Ellis and J. I. Kapusta, Phys. Rev. C 64, 035202 (2001) 

[12] A.Dobado, A.Gomez Nicola, F.J.Llanes-Estrada and J.R.Pelaez, arXiv:hep-ph/0206238. 

[13] T. N. Truong, Phys. Rev. Lett. 67, 2260 (1991); A. Dobado, M. J. Hcrrero and T. N. Truong, Phys. 
Lett. B 235, 134 (1990); A. Dobado and J. R. Pclacz, Phys. Rev. D 56, 3057 (1997). 

[14] M.Le Bellac, Thermal Field Theory (Cambridge University Press, 1996) and references therein. 

[15] R.Baier and A.Niegawa, Phys.Rev. D 49 (1994) 4107. 

[16] R.Kobcs, Phys.Rev.D 42 (1990) 562; 43 (1991) 1269. 

[17] A.Niegawa Phys.Lett. B 247 (1990) 351: Phys.Rev. D 57 (1998) 1379. N.Ashida, H.Nakkagawa, 
A.Niegawa and HYokota, Phys.Rev. D 45 (1992) 2066; Ann.Phys. (NY.) 215 (1992) 315; 230 
(1994) 161 (E). 

[18] H. A. Wcldon, Phys. Rev. D 28, 2007 (1983); Annals Phys. 214, 152 (1992). 

[19] J.F.Donoghue, C.Ramirez and G.Valencia, Phys. Rev. D39, 1947 (1989); G. Ecker, J. Gasser, 
A. Pich and E. dc Rafael, Nucl. Phys. B 321 311 (1989). 

[20] C. A. Dominguez, M. Loewe and J. C. Rojas, Z. Phys. C 59, 63 (1993); R. D. Pisarski, Phys. Rev. 
D 52, 3773 (1995). 



11 



